TSCHiifidAL MEKOEAiTDUMS 
NATIONAL ADVISOSY (JoiiMlilTES S-OH Afi P.O ST AUT I C S 



ITo. 972 



USE OJ CHARTS EOE PLOW DISCHAEaS CALCULATIONS 

3y 0. LMtz 



L"af t fahrtf orsciiung 
Yol. 17, Ko. 8, Octo"ber 26, 1940 
Verlag von R. OldentoTirg , M-anchen und Berlin 



Trashing ton 
March 1941 



mil 




Hill 



3 1176 01440 7200 



KATIOVjlL ASTISOBT commutes TOB AEBOirAUTICS 



. . . lECHiriCAL UEMORANDUH KO. 972 



USX 07 CEA.BTS 70 B TL07 DISCEjLBaE CALCULATIONS* 

By 0. LtitE 



SUMIUBY 



Yarloua problems In connection vltta engine design 
Involve f lov-dlscbarge calculations which are rendered 
difficult hoth on account of the large numher of external 
varlahlea that enter Into the computation - I.e.. changes 
In discharge area during the process, change In volume of 
the cylinder, pressure, etc., and changes In the thermal 
constants themselves of the flow medium. a fairly accu- 
rate solution that does not Involve an excessive amount 
of lahor can he obtained only through the extensive use 
of 1-B tables. In the present report, a- solution la of- 
fered In the form of a different method making use of the 
I-S tables of Lutz and Wolf. 



1. INTBODUCTIOS 



Those flow processes of most Interest In connection 
with engine design occur mainly In the supercritical 
region at temperatures for which the specific heats, and 
also the adlabatlc coefficient, must no longer be con- 
sidered as constant. Tor these conditions l-s charts are 
employed. The simplest and most convenient for this pur- 
pose are the I.-S. charts of 0. Lut s and T. Wolf. On ac- 
count of the large number of external variables Involved - 
discharge area, cylinder volume, back pressure, etc. - a 
straightforward numerical solution of the problem Is Im- 
possible and attempts In this direction lead only to com- 
plications (reference 1). 

A different method Is -here developed which has the 
advantage of clearness and which takes Into account all 
the variables that enter the problem. 



Die graphlsche Irmlttlung von Ausstromvorg«ngen. Luft- 
f ahrtf orsehung, vol. 17, no. 8, Oct. 26. 1940. pp. 
332-335. 
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2. DEBITATiaiir 07 !rH3 BISLATIOlirs US15I) 



The discharge 1b assumed to talre place from a ^working 
cylinder (fig* 1) whose ▼olume at any instant is deter- 
mined "by the piston area -f^.' and the displacement h of 
the piston from the cylinder head. With carved cylinder 
head or piston head, the value of h is given hy the con- 
dition that the ratio of h^^^ to the piston stroke h^ 

la equal to the ratio of e to e - 1 where c is the 
compression ratio. Ine dlscnarge area (taking Into ac- 
count also constrictions in the Jet) will he denoted hy 
f, the discharge velocity at the minimum Jet cross sec- 
tion by Wq. Magnitudes characterising .the state of the 
gas in the cylinder will have tne suhscrlpt 1. 

For Wj^, we have. 

Wa = 91.5yi:lJl-i: (1) 

where Z. is the Internal heat or enthalpy in the dis- 
charge area and U the gram molecular weight of the dis- 
charge medium. In the Z,-S tables referred tOi(page 7). 
the method is indicated of obtaining these values for com- 
bustion gases whose thermal properties deviate from those 
of air. The Internal heat I, in the discharge section 
will be equal to that in the external region only for the 
subcritlcal flow condition; for the supercritical condi- 
tion it must be separately determined. 

Assuming the process is adiabaticiwe have for the 
change in state in the cylinder 

^Pl _ £pi dC^i ■ dh 

where C. denotes the gram molecular heat, the instan- 

taneous weight of the gas and dh/h the relative change in 
volume. The time rates of change are connected by 



-Lk = SeI Pi %)i Pi 
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ZTow dh/dt Is the InstantaneoxiB piston velocity o and 
dOj^/dt may te replaced ty - ^fw^ ap t)iat 



dt b V^fj. Yi • > 



which, may also he written 



dp^ S.pi P 



f 7 



dt .X i 

where T le the gram molecular Tolume. 

It la convenient in engine comp'itat ionfl to refer to 
the crank angle qp. Since the time element dt and the 
crank-angle element dq? are connected hy the relation 

dV 

dt - -— (n = rpm) 
on 

we have finally^ paaalng; to finite differencea Apj^ 

The known external variables are the revolxitlona per 
minute n, the distance of the platon from the cylinder 
head h, the ratio of croaa aectlona t/f-^, and the In- 

atantaneoxLB platon velocity c; the Initial atate p^, 

7.4 muat furthermore he given. The ratio Ki = /c 

1 "^1 ~^1 

of the apeclflc heata and the ratio Z^/l. of the gran 

molecular volumea for aupererltleal flowa muat be de- 
termined by further eonalderatlona* 

In figure 3ithe val,ues of k are plotted as func- 
tions of the temperature with the factor p of the 1,-S, 
tables (page 5) as parameter. It may be pointed out in 
this connection that the valuea of k can alao be 
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dlrectl7 determined frofii the 1-1 tables If net too great 
accurac7 Is required. lor the constant-pressnre cnrves 
of the IrS, tables we have, as is known, 

ds = c i?. 



and for the constant- volume curves 




Again passing to finite differences, we obtain for the re- 
quired ratio of the specific heats 




Accordingly, to find Cp/2y, the increments AS^ and 
between the isotropea. isobars and Isochores (con- 
stant-volume curves) corresponding to a small increment 
AT from the state a, are read off on the L-S, chart as 
Indicated in figure 3. We then have 




(4a) 



To reduce the error due -to the curvature of the curves, 
it is best to carry out the process for positive and 
negative temperature differences and average the results. 

The exact determination of the critical pressure 
ratio from which the discharge velocity for supercritical 
flows and the corresponding molecular volume are obtained 
is not possible in explicit, form. Assuming that, for 
supercritical flows, the weight passing through the dls- 

charge area f =— remains a maximum, we must have 
— s 
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or using (1) 

' 2(ii - L3) 



~ - =..P. .. . . (5) 



"So 



V dl„ = C- dT. and « — for an adla'batlc process 



= - — S. BO that we have finall7 



AH Tb 

This equation, which as yet has not heen simplified by 
neglecting terms, can he solved only graphically. The 
temperature ratio Tg/Tj. was determined from the I-S chart, 
nuaher 1, 1* heing determined from the state I... T_. 

Here. too. an approximate method can he indicated. 
Equation (6) can he simplified by assuming the specific 
heat for the expansion process from the inside of the 
cylinder to the opening to be constant and equal to the 
▼alue Cp In the opening. 7e then have 



Ts 2CV3 



(6a) 



The right-hand side is determined as in figure 3 except 
that the segment IZ is taken only up to the midpoints be- 
tween the lines of constant volume and the isobars, as 
shown in figure 4. We then have 



- Tb Z 

Tl'll <6^) 

Thus Tg and therefore the state in the I.- 3, chart can be 
det ermlned. 
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Tor the ease of flovs below the erltloal oondltlon, 
the final state Is that of the external region and hence 
known. 

All the magnittideB In formula (3a) are now Icnown and 
the discharge can he cooiputed step hy step without too 
much lahor. 



The discharge process from a cylinder of 1.64 liters 
volume will he determined, the initial state being given 

Pi - °»0 atm, *i = 1400°; the valve opening time Is 
shown In figure 5. According to equation (3a) for the de- 
termination of the pressure variation, there are first re- 
quired to he known the values h, c and f/fj^ associated 

with the engine dimenBlons. The variation of these values 
with the crank angle is shown in figure 5. Since f /f ]^ 

denotes the ratio of the minimum Jet cross section (not 
the valve cross section f^) to the piston area, it is 

necessary hy a separate test to determine the contraction 
coefficient of the valve under consideration. In fig- 

ure 6, the values of \k for a few valve arrangements are 
plotted.* Corresponding to tnese data, curve 4 has been 
chosen as approximately applying to our case and f/f]^ 

determined in figure 5. 

It was found convenient to use steps of 10° crank 
angle for the supercritical discharge region and 5° for 
the suberltlcal. 

For the first step, the initial values are -p^ = 6 atm. 



!*1 = 1400°, hggo = 0.141 m, 05^0 = 15.63 m/s and ^ -f- 



= 0.00083. The ratio of the specific heats at the 



According to tests conducted by E. Hummel at the Instltut 
fur Uotorenf orschung der Luf tf ahrtfor schungsanstalt 
Hermann aorung. In applying statically measured values 
'to dynamic processes, it should be remembered that the 
value correction factors may vary as a result of pos- 
sible pulsations in discharging. Account will be taken 
of this fact in a later report. 



3. 



ILLUSTRLTITZ EXAICPLE 





* 
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Initial testperature Tj^ = 1673° E 1b found from figure 2 

to have the value 1.232 (where, for the comhuetlon gas, a 
stoichiometric gasollne-alr mixture of tahle factor p 
= 1.50 Is to he assumed, see 1,-S charts). Trom figure 2, 
the temperature ratio T^/Tj^ 1b obtained aa 0.882 with 

the aid of which the molecular volume ratio Vj^/Z, is 

found from the 1^-8, charts to have the value 0.612. The 
discharge veloclt7 w^, is finally obtained as 743 m/s. 

We thus have for the first step: 

Ap = - 1.262 -^-5 (0.00083X0.612X743 + 15.63) X 10 

6X2700X0.141 

= - 0.530 atm 

This pressure drop is large so that the first step (for 
which the Initial values of 6.0 atm and 1400° were used 
and not the average values which should have heen suhst i- 
tuted for more accurate computation) must be corrected. 
The computation is therefore repented with pj^ = 5.74 atm 

and = 1384° (from tne I.-5 charts) and there is ob- 

tained 

5 7 40 

Ap = - 1.262 2 70*0X0. 1 41 (0.00083X0.620X783 + 15.63) X 10 

- - 0.508 atm 

There is t£.en obtained for the next step the initial value 
Pl50° ~ ^•^9'^ B.tta. and from the I-S charts ^j^sqO = 1369°. 

These values are indicated in figure 7. It is now possible 
with sufficient accuracy to estimate the mean computation 
value Pi460 for the second step. Generally no correc- 
tion will be required. 

The accuracy of the computation procedure depends on 

a) the sise of the computation steps, 

b) the accuracy with which the values . v**-^* 

5.V1 ^ I 

and Wq (see equation (3a)) can be estimated 
or corrected before each half step, 
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c) the accuracy in reading off the I'i. charts In 

determining values and w^^. It la "better 

in place of T to obtain the corresponding 
▼aluBB p and ^ which can he more accurate- 
ly read off and then determine 7 from the 
gas equation. 

Figure 7 shows the results of the computation. On 
figure 8 is shown finally the weight of gas 0* dis- 
charged per unit of time. The curve has a maximum short- 
ly after reaching inner dead center and toward the end of 
the process has a slight irregularity which is due to the 
large piston velocity. 



Translation by S. Heiss, 
National Advisory Committee 
for Aeronautics. 



1. Schule, 7.: Technischo Thermodynamik. Ed. 2, S. 390 ■ 
Berlin, 1923. 
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figs. 1,2,3, 4 
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Figure 2.- Values of k ani as functions of the teinperature and the 



chart factor 0 . 



Ssconst. vsconBt. 
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i=con8t. 



Vi.gMxe 3.~ Determination ofK. 



Tigore 1.- Sketch indicating 
notation. 




Jigure 4.- Detennination of _JI 

Ti 
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Iig».6.6,7,8 




Vlgurc 5.- Platon dlsplacoaent h., 

piston Velocity c, 
valT* cross section, fy .contrac-' 
tion coefflciantilr and ratio f/fjc 
of the discharge area to the 
piston area. 



Figure 6.~ Contraction coeffi- 
cient ^ for various 
valve arrangements. 
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Figure 8.- Discharge weight 0* o.ob 
and valve lift hy o.oe 
as a function of the crank o.a^ 
angle position. 
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